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ABSTRACT 

We study the effect of the supersonic baryon-CDM flow, which has recently been shown to have a 
large effect on structure formation during the dark ages 10 < z < 1000, on the abundance of luminous, 
low-mass satellite galaxies around galaxies like the Milky Way. As the supersonic baryon-CDM flow 
significantly suppresses both the number of halos formed and the amount of baryons accreted onto 
such halos of masses 10 6 < M^q/Mq < 10 8 at z > 10, a large effect results on the stellar luminosity 
function before reionization. As halos of these masses are believed to have very little star formation 
after reionization due to the effects of photo-heating by the ultraviolet background, this effect persists 
to the present day. We calculate that the number of low-mass 10 6 < M^iq/Mq < 10 8 halos that 
host luminous satellite galaxies today is typically suppressed by 50 percent, with values ranging up 
to 90 percent in regions where the initial supersonic velocity is high. We show that this previously- 
ignored cosmological effect resolves most of the tension between the observed and predicted number of 
low-mass satellites in the Milky Way, obviating the need for any other mass-dependent star-formation 
suppression before reionization. 

Subject headings: cosmology: theory — early universe — galaxies: formation — Galaxy: halo - 
galaxies: statistics — Galaxy: structure 



1. INTRODUCTION 

A robust prediction of the Cold Dark Matter (CDM) 
cosmological paradigm of hierarchical clustering is that 
the halos of galaxies like the Milky Way should contain 
hundreds of satellite subhalos that could be expected 
to host observabl e galaxies (Kauffmann et al.l 119931 : 
iKlvpin et all 119991 : iMoore et al.l 119991 ). T he success of 
the A CDM model from the horizon-scale (iKomatsu et al 



discovered many more smaller dSphs in the Milky Way's 

halo (e.g., ^^^n^^_^^_E ^^^ ^^_^M 

iZucker et all 2006; Belokuro vet al.l 120071 : llrwin et all 



2011 ) to the small scale of the Lyman-a forest (jViel et al 



2008ft suggests that this "missing satellites problem" is 



most likely the result of the details of galaxy-formation 
physics in low-mass halos rather than a manifestation of 
a deviation from the standard framework on small scales. 

Potential solutions to the missing satellites problem 
have so far come in two flavors. One approach aims 
to reduce the intrinsic small-scale power in the CDM 
framework by positing alternatives to standard infla- 
tionary mechanis ms for producing the initial pertur- 
bation spectrum (jKamionkowski fc Lld dle 2000). Sim- 
ilarly, warm dark matter models naturally suppress 
the number of small, boun d structures in the Universe 
([Polisenskv fc Ricottil 1201 J ). The alternative to mod- 
ifying the cosmological framework is to invoke astro- 
physical explanations for suppressing star formation in 
low- mass galaxies. In particular, suppression of star 
formation by photo-heating by the ultraviolet (UV) 
background after reionization naturally leads to a star- 
formation-efficien cy cut-off at approximately the cor - 
rect mass scale (B ullock et all 120001 : iSomervilld [2002ft . 
This can explain the number counts of "classical dwarf 
spheroidals" (dSphs) around the Milky Way (My < —6; 
IKodosov etaLl 12009ft . 



The Sloan Digital Sky Survey (SDSS: lYork et al.l l2C 



Hubble fellow 



120071 : IWalsh et al.ll2007ft . To explain the number counts 
of these new discoveries within ACDM, it is neces- 
sary to invoke additional suppression of star forma- 
tion in the lowest m ass halos before reionization (e.g., 
iKoposov et al.l 12009). as star formation in these ha- 
los mostly ends at reionization due to photoioniza- 
tion. Specifically, the stellar mass in low-mass halos 
(M < 10 8 M^) can only be < 10 -4 times the halo mass 
(jMadau et al.l l2008h . an order of magnitude less than 
the stellar- mass fraction in higher mass halos. This is 
only w 0.05 % of the universal baryon fraction, while 
typically we would expect star formation efficiencies 
of at least a few percent, even when star formation 
only proceeds through iif 2 cooling ([Bovill fc Ric otti 2009; 
Salvador i fc Ferra ra 2009) and it is unclear whether ra- 



diative feedback from the first generation of stars sup- 
presses pre- reionization star formation in low-mass ha- 
los or not (lHaiman et all 119961: [Omuk ai & Ni shjl H999t 



Hai man et al. 2000; Ricotti et al. 2002a : b; Whalen et al. 
2Q08D. 

In this paper we determine the influence of the 
previously-ignored effect of the relative velocity between 
baryons and dark matter at recombination on the stel- 
lar content of the smal lest-mass galaxies. As was re - 
cently pointed out by iTseliakhovich fc Hiratal (|2010ft . 
the typical baryon-CDM relative velocity of 30 km s _1 
at recombination is supersonic as the baryon sound 
speed at kinetic decoupling drops to 5 km s _1 . While 
the relative velocity decays as ex a -1 , where a is the 
scale factor, the supersonic baryon flow has a pro- 
found effect on the formation of the first structures 



(ITseliakhovich et al.l 120111 iStacv et alJl20lH iMaio et al 
12011b iGreif et al.l 120111 : iNaoz et a l. 2012; Visbal et al 
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120 121 : IQ'Learv fc McQuinril2012l ). which through the ef- 
fects of hierarchical clustering m ight persist today (e.g., 
in the baryon acoustic feature, iDalal et al.ll2010l ). and 
it could also give rise to CMB B-modes if the effect 
persi sts during the epoch of reionization (Fe rraro et al.l 
l2Q12f ). We show that the supersonic baryon flow has a 
large effect on both the number and the baryon content 
of halos with masses between the H2 cooling limit (~ 
1O 6 M ; iTegmark et al.l 119971) and the p hoto-ionization 
scale (« 1O 7 M at z = 11. IGnedinll2QQQl ). This typically 
reduces the number of satellites by 50 up to 90 percent 
at My « -1. 

The outline of this paper is as follows. In §[2] we review 
the supersonic baryon-CDM relative velocity effect. We 
then compute the pre-reionization mass and luminosity 
functions in § O We discuss the low-redshift behavior of 
the satellite luminosity function in § [U Our conclusions 
are in §0 We assume cos mological parameters matching 
the WMAP 7-year data (jKomatsu et all 120111 ): ft 6>0 = 
0.0456, ft c>0 = 0.227, z eq = 3232, H = 70 .4 km s" 1 
Mpc" 1 , A 2 ((k = 0.002 Mpc -1 ) = 2.44 x 10~ 9 , and n s = 
0.963. 

2. LINEAR EVOLUTION OF DENSITY PERTURBATIONS 
IN THE PRESENCE OF A SUPERSONIC FLOW 



Aft er kinetic decoupling (z ~ 1020; lEisenstein fc Hul 
119981 ) the baryonic sound speed drops to 5 km s _1 , while 
the baryons move relative to the CDM with a typical 
velocity of 30 km s -1 . As originally pointed out by 
Tseliakh ovich fc Hiratal ([2010), this means that second- 
order terms such as v • V5 and v • Vv, which are typ- 
ically neglected to first order, become as large as first- 
order terms in the continuity and Navier-Stokes equa- 
tions that describe the evolution of inhomogeneities in 
the baryons and CDM after recomb i nation . We follow 
the treatment of iTseliakhovich et al.l (|2011l ). who shows 
that the supersonic relative flow Vb c is homogeneous over 
a few comoving Mpc, with Vb c drawn from a Gaussian 
with a variance per axis of <J 2 c /3 with &bc = 30 km s _1 at 
kinetic decoupling. The evolution equations for the rela- 
tive CDM (S c ) and the baryon (£5) perturbations can be 
written in the baryon rest frame in a ^ comoving Mpc 
patch with a supersonic flow V5 C as 



— = -v 6c .kd c - C 
at a 



6,0 
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where we write the baryonic sound speed as c 2 s = 
kBTb(5b + $T h ) I ' iirriH^b-, is the velocity divergence, 
^m,o = ^6,0 + ^c,o? M = 1-22 is the mean molecular 
weight including a helium mass fraction of 0.24, run 
is the mass of the hydrogen atom, and X5 and Sr b are 



the mean baryon temperature and its relative fluctu- 
ation, respectively. The evo lution of X5 is given by 
(jTseliakhovich fc Hiratal I2010D 



^CMB,0 



1 



q/(l/119) 
l + (l/115/a) 3 / 2 



(2) 



where Xcmb,o = 2.726 K. The temperature fluc- 
tuation evolves accor ding to (jBarkana fc Loebl 120051 : 
iNaoz fc Barkanall2005f ) 



9S Tb 
dt 



--Vh 

3 6 



x e (t) 1 T 7 ^ 
t-y a 4 Tb 



(3) 



where x e (t) is the free electron fraction out of the to- 
tal number density of gas particles, t 1 = 8.55 x 10 -13 
yr _1 , T 7 is the mean photon temperature, and we 
ha ve neglected the photon inhomogeneities compared 
to INaoz fc Barkanal ([20051 ). We ob tain the free elec- 
tron fraction x e (t) using RECFAST (jSeager et al.lll999l 
[2003). 

We solve the complex system of Equations (p~H3]) by 
taking initial conditions for the matter inhom ogeneities 
and velocities from CAMB0 ( Lewis et al.ll2000l ) at kinetic 
decoupling (z = 1020) and setting up the initial inho- 
mogeneities in the b aryon temperature follow ing the ap- 
proximation done in INaoz fc Barkanal (|2005l ) by requir- 
ing 8ST b /dt = 05^ /dt, where 5^ is the relative pho- 
ton temperature perturbation. This approximation is 
justified because of tight thermal coupling between the 
baryons and the photons and it has been shown to af- 
fect the power spectra only by a fraction of a percent at 
z = 1020. The relative velocity in Equation (pQ) decays 
as Vfo c oc a -1 . 

3. LUMINOSITY FUNCTION BEFORE REIONIZATION 



We follow the procedure of ITseliakh ovich et al.l (|2011f ) 
to calculate the halo mass function and the baryonic con- 
tent of low-mass halos before reionization in the pres- 
ence of a supersonic baryon-CDM flow. We calculate the 
halo mass function using th e Extended Press- Schechter 
forma lism (jBond et a,lJll991b lBowerlll99H : lLacev fc Coll 
fl99l 



dN 



dMi 



halo 



■(M ha i ;v bc ) 



Po 



alo 



dS 



dM, 



halo 



f(S c (z),S), (4) 



where po is the mean matter density in the universe and 
S is the usual variance given by 



S(M;v bc ) = / d\nkA 2 m (k;v bc )\W(k;R)\' 2 



(5) 



Here A^(/c;^6 C ) is the matter power spectrum com- 
puted using Equations (p~H3]). The transfer functions 
computed using Equations ([THH]) depend on the angle 
between k and V5 C . In what follows, we average all so- 
lutions for \Sb\ and \S C \ over this angle and give results 
for Vbc = |vfc c |. W(k]R) is the tophat window function 
and the initial comoving radius R corresponds to a halo 

2 |http : //camb . info/ 1 . 
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Fig. 1. — Halo mass function at z = 11 with a relative baryon- 
CDM velocity of Vb c = and Vb c / 0. 1 at, c corresponds to approxi- 
mately 30 km s _1 at decoupling [z = 1020). The globally-averaged 
effect obtained by averaging over the Gaussian probability distri- 
bution function of v^ c is shown as the gray curve. 

of m ass Mhaio = 47ri? 3 po/3. We use the functional form 
from lSheth fc Tormenl (fl999h for f(5 c (z),S): 



f(5 c {z) 1 S) = A^J^- 



-av 2 /2 



(6) 



where v 6 n (z)/y/S \ a = 0.75, q = 0.3, A = 0.322 
(jSheth fe Tormenl I2002D . and S c (z) = 1.67 is the critical 
density of collapse at z ~ 10. 

Constraints from the Cosmic Microwave Background 
on the epoch of reionization are consistent with a red- 
shift of reioniza tion of z = 11 (K omatsu et al.l 120111 : 
IZahn et"aTll2011f ) (defined as the redshift at which reion- 
ization would begin if the universe was reionized instan- 
taneously). The halo mass function at z = 11 cal- 
culated using the procedure given above is shown in 
FIG.[TJ This figure shows that the effect of the supersonic 
baryon-CDM velocity typically suppresses the number 
of dark matter halos by 10 to 20 percent in the range 
10 5 < M halo /M^ < 10 7 at z = 11. TV-body simulations 
bv lNaoz et al.l (|2012[ ) show a similar level of suppression 
in the mass function at z = 11. 

To estimate the stellar mass in low- mass subhalos at 
reioni zatio n, we use the res ults from iNaoz fc Barkanal 
(|2007l ) and INaoz et al.l (|2009l ). These studies have shown 
that the gas fraction in halos at high-redshift is sup- 
pressed with respect to the universal baryon fraction 
by the combined effect of the remaining suppression af- 
ter recombinati on in the baryon density perturbations 
on sm all scales ([Naoz fc Ba rkana 2007; iBarkana fc Loebl 
|2011[ ) and the baryonic pressure, which gives rise to a 



Fig. 2. — Luminosity function at z = 11 using the luminosity 
that the halos would have today. This is the total effect on the 
luminosity function from the combination of the suppression of 
the halo mass function and the reduction in the gas fraction due 
to Vb c / 0. The average case is shown in the gray curve. The 
separate contributions to the suppression for v^ c / from the halo 
mass function and gas fraction are shown in FIG. [3] The arrow in 
the bottom panel shows how far each curve moves when lowering 
the star-formation efficiency to f s = 0.005. 

redshift- dependent ' filtering' scale jjGnedin fc Y. 



Gnedin 2000; Naoz & Barkana 2007). This filtering scale 
corresponds to the length scales below which the baryon- 
to-total-matter fluctuation drops substantially below its 
large-scale value, and simulations find that it charac- 
terizes the minimum halo mass that can ret ain its gas 
(|Naoz et alJl20Q9HO ? Leary fc McQuinnll2012D . To deter- 
mine the filtering scale kp as a function of the supersonic 
flow velocity V5 C we fit the functional form 



1^6 



\St, 



n 1 + r L ss 



(7) 



where |£ tot | = (^5,0 + ^c,o |£ c |)/^m,o, ^lss is deter- 
mined from the behavior of |£b|/|#tot| on larger scales, 
1 < k Mpc < 10, and kp and n are subsequently fit at 
k > 1 Mpc -1 . 

Using the filtering scale fc p, the filtering mass is defined 
as (|Naoz fc Barkanal 1 2 00 fh 



M F (M halo ;v bc ) 



47T 



-Po 



7T 



(8) 



Using the filtering mass, we can estimate t he gas frac- 
tion at redshifts z > 11 as (INaoz et al.ll2009l) 



f gas (M;v bc ) = f b0 



l + (2 
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Fig. 3. — Separate contributions to the difference in the z = 11 
luminosity function in FIG. [2] between Vb c = and Vb c ^ from the 
suppression in the halo mass function in FIG.^and the suppression 
in the gas fraction. 

where /^o is the gas fraction in the high-mass limit. 
fb,o is typically lower than the universal baryon frac- 
tion because density perturbations in the baryons at 
high redshift remain suppressed due to the lingering ef- 
fect of the coupling b etween baryons and pho tons be- 
fore recombination. In Bark ana fc Lo"ebl ([2011] ), the au- 
thors show that fb.o ~ (1 + 3.2rLss)^fr.o/^m.o- Following 
Tseliakh ovich et all (|2Q11[ ), we use a = 0.7, even though 
this value of a was calculated for halos at z ~ 20. As 
our main objective is to show the difference between the 
Vbc — and Vb c ^ cases, the exact functional form in 
Equation (j9]) does not matter greatly. 

We assume that a fraction f s of the gas in a halo — 
independent of halo mass — turns into stars, such that 
the stellar mass is given by 

M s (M ha i ; v bc ) = fs / g as(M ha i ; v hc ) M hsXo . (10) 

We convert this stellar mass into a luminosity by assum- 
ing that this stellar mass is turned into stars with very 
low metallicity Z = Z /2OO. Today, such a population 
would shin e with an absolute magnitude My = 6.7 per 
solar mass (jMadau et al.l 2008: B ruzual fc Charlotll2003l ) 
and it is this luminosity that we use to show the resulting 
luminosity function. 

FIG. [2] shows the luminosity function at z — 11 using 
the luminosity that the halos would have today. That is, 
it is the luminosity function we would observe today if 
the low-mass halos were frozen in their pre-reionization 
state and there was no evolution in the number of halos. 
We will study their evolution in § 3J 

FIG. [2] assumes a star formation efficiency of f s = 0.03, 
which gives a total-to-stellar mass at the high-mass end 
of 10 -2 - 8 . The arrow in the bottom panel of FIG. [2] 
shows the effect of lowering the star-formation efficiency 
to f s = 0.005. In this case, halos of a given luminosity 
come from higher total-mass halos, where the effect of the 
supersonic baryon-CDM velocity is smaller, thus shifting 
the effect toward smaller luminosities. 

FIG. [3] shows the relative contribution to the total ef- 
fect in FIG. [2] from the suppression of the mass function 
in the presence of the supersonic baryon-CDM flow (see 



FIG. [I]) and the suppression of the baryon fractions. It is 
interesting to note that the suppression in the accretion 
of baryons onto dark matter halos is the main driver of 
t he suppression of low- l umino sity halos. 

iMcQuinn fc O'Learvl ()2012h find that shock heating 
can raise the baryon temperature by approximately 
10 percent compared to the evolution that we assume. 
If we approximate the evolution in that they find by 
raising T5 by 10 percent at z < 20, we find a negligible in- 
fluence on the matter power spectrum and filtering mass 
at z = 11, such that the conclusions of this section are 
unaffected by shock heating. 

4. PRESENT-DAY SATELLITE LUMINOSITY FUNCTION 

In this section we quantify how the suppression at the 
faint end of the luminosity function before reionization 
affects the present-day subhalo luminosity function for 
a Milky- Way sized halo. We estimate the satellite lumi- 
nosity function by running merger-tree si mulations using 
the Ex tend ed Press- Schechter formalism of Lac ev fc Cold 



([19931 ) and iCole et al.l (|2000l) fo r a parent halo of mass 
M = 10 12 M© (jXue et aljl2008h . We create 100 merger 
trees for each of Vb c = 0,1,2, and 3 &bc by computing 
the z = S(M; Vb c ) as before and using the linear over- 
density for collapse S c = 1.67, extrapolated to z = 
using the growth factor (Ek e et all 119961 : iCarroll et al.l 
Il992f ). We resolve subhalos down to a resolution mass 
of M res = 10 6 M© and run the simulation until z = 11. 
As we are only interested in low-mass subhalos, we only 
track the mass evolution of subhalos, without resolving 
it into sub-subhalos. For each subhalo, we record the 
pre-reionization mass M(z = 11), defined as the mass of 
the subhalo at z = 11, and the mass at the time of the 
merger with the parent halo M(z = z acc ). 

We do not follow the dynamical and mass evolution 
of the subhalos after they merge with the parent halo. 
This approximation assumes that low-mass (< 1O 8 M ) 
satellites do not get tidally disrupted and do not lose a 
significant amount of stellar mass after merging with the 
parent halo. High-resolution iV-body simulations find 
that only a few percent of M > 1O 7 M0 subhalos are 
tidally destroyed between z = 1 and z = and that sub- 
halos with M < 10 8 M ^ retain most of their total mass 
([Diemand et al.l 12007). Even when the outer parts of 
the subhalos are tidally stripped, the stars and the inner 
part of the subhalo are stripped only at the last stage 
of th e tidal disruption of the subhalo (P eharrubia et al.l 
l2008f h such that most of the stellar mass is retained, 
even if a significant part of the dark matter halo is 
tidally stripped. Therefore, to a good approximation, 
the stellar-mass function of low-mass satellites should not 
be strongly affected by the effects of tidal stripping and 
disruption. 

We compute the luminosity of each subhalo by us- 
ing the prescription of Equation ([TO]) applied to the 
pre-reionization mass M(z = 11) to calculate the pre- 
reionization stellar-mass and turning this into a lumi- 
nosity again using My = 6.7 per solar mass. We cal- 
culate the post-reionization stellar mass for each sub- 
halo by using a star-formation efficiency t hat takes into 
account suppression by photoionizati on ([Gnedinl 120001 : 
iHoeft et al. 2006: Qk amoto et al. 2008) and that assumes 
that the effect of the supersonic baryon-CDM velocity 
has no influence on star formation after reionization: 
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Fig. 4. — Satellite luminosity function at z — for a Milky- Way 
size halo (M = 10 12 Mq). The curves show the average luminosity 
function of 100 merger trees for each value of Vb c . The bottom 
panel shows the fractional effect with respect to Vb c = 0. 



M s (z < 11) = /„ 



(M(z = z acc ) - M(z = 11)) 
(l + 0.26 (V ait /V chc (z acc )Y 



> (11) 



where = 10 3 / 6.25 and V cr it = 35 km s 1 (as in 



IKoposov et al.ll2009j ). We calcul ate V c \rc(z H , cc ) by usin g 
the virial radius of Equation (1) in lKoposov et al.l (|2QQ9l ). 
The post-reionization contribution to the total stellar lu- 
minosity of a satellite is then calculated by assuming a 
solar mass-to-light ratio (Ma rtin et al.ll2008l ). 

The luminosity function of satellites of a Milky- Way 
type halo computed in this way is shown in FIG. |U A 
comparison between this figure and FIG. [2] confirms that 
the pre-reionization effect on the suppression of the lu- 
minosity function at z = 11 translates into almost the 
same suppression in the satellite luminosity function at 
z = 0. Additional suppression of the luminosity function 
results from the fact that in a region of a highly super- 
sonic baryon-CDM flow, structure formation is delayed, 
such that a larger fraction of the present-day mass of a 
halo was accreted after the epoch of reionization, where it 
is affected by the suppressed gas accretion and cooling in 
the presence of the photoionizing background after reion- 
ization. This additional effect also explains the slight in- 
crease in the number of brighter satellites in FIG. SJ At 
the bright end — where pre-reionization suppression due 
to the supersonic baryon-CDM flow is small — a larger 
fraction of a satellite's mass is accreted after reioiniza- 
tion where we have assumed a mass-to-light ratio appro- 
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Fig. 5. — Luminosity function of satellite galaxies of the Milky 
Way observable by t h e SPS S. The diamond data points are taken 
from Koposov et al. (2009). The error bars on the Vb c = la^c 
model curve show the 68 percent spread in the merger trees at 
luminosities where the predictions with different v^c differ; this 
spread is similar for all four model curves. The error bars are highly 
correlated; for example, the correlation between the My = —4 and 
My = —2 error bars is 0.75. For this figure we assume / s = 0.01. 

priate for a metal-poor population rather than that of 
an extremely metal-poor population, which we assumed 
for pre-reionization star formation. This discontinuity, 
which is an artifact of the simplicity of our assumptions, 
stretches the ^ curves toward the brighter end, 
leading to an increased number of bright satellites. 

To compare the predicted satellite luminosity function 
to the observed l uminosity function o f satellites of the 
Milky Way (e.g., IKoposov et all 120081 ). we assume that 
the spatial distribution of the satellites follows a Navarro- 
Frenk- White profile with a scale radiu s of 30 kpc; in 
accor dance with numerical simulations ([Diemand et al.l 
we anti-bias this distribution radially by multi- 
plying by the Galactocentric radius. We then calculate 
the observational fraction by using the si mple model for 
the S PSS satellite selection function of IKoposov et al.l 
(|2009l ). where a satellite of absolute magnitude My can 
be detected out to a distance of Anax = lO 11 " 228 ^ 
kpc from t he Sun (assumed to be 8 kpc from the Galac- 
tic center, iBovv et ai]|2009f h integrating out to a virial 
radius of 260 kpc. We multiply this selection fraction 
by 0.194 to account for the partial sky coverage of the 
SDSS. All satellites down to My & —6 could be observed 
throughout the virial volume by the SDSS. 

The luminosity function of satellites detectable by the 
SDSS thus computed is shown in FIG. El Even though, 
for this figure we have assumed a star-formation effi- 
ciency f s of only 0.01, it is clear that the supersonic 
baryon-CDM velocity has a large effect on the faint end 
of this luminosity function that significantly lowers the 
predicted fraction of satellites at My > —6. Thus, the 
inefficient gas accretion at high redshift induced by a 
large supersonic flow has a large and lasting effect on 
the abundance of luminous satellites for a galaxy like the 
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Milky Way. 

We do not show the upper limit of dN/dMy < 0.1 at 
My = from lKoposov et all (|2009D in FIG. El as the pre- 
dictions at My = strongly depend on the exact form 
of the SDSS selection function. At My = 0, the maxi- 
mum distance out to which a satellite can be detected by 
the SDSS is approximately 12 kpc. As we do not follow 
the dynamical evolution of the satellites in the parent 
halo, our simulations do not capture the fact that those 
satellites that come within 20 kpc of the Galactic cen- 
ter have a much larger chance of being tidally disrupted, 
such that our predictions at My = over-estimate the 
number of observable satellites. 

5. CONCLUSION 

As pointed out by Tselia khovich fc Hiratal (|2010f h at 
recombination, baryons move with a typical supersonic 
velocity of 30 km s _1 with respect to the dark mat- 
ter. While this relative baryon-CDM velocity decays as 
oc a -1 , it has a large effect on the formation of the first 
structures in the Universe. We have investigated the ef- 
fect on the present-day abundance of luminous low-mass 
satellite galaxies (< 1O 8 M ) of a Milky- Way size galaxy. 
One would expect the supersonic baryon-CDM veloc- 
ity to affect the abundance and luminosity of low-mass 
satellites, as these satellites must form many of their 
stars before reionization because photo-heating by the 
UV background after reionization suppresses star forma- 
tion in halos with circular velocities < 30 km s _1 . The 
effect of the supersonic baryon-CDM flow should be a 
standard cosmological effect that needs to be taken into 
account when predicting the number of luminous satel- 
lite galaxies — the only unknown being the initial value of 
the relative velocity in the patch of the Universe that is 
being considered. Additionally, we are motivated by the 
apparent discrepancy between the observed number of 
dwarf satellite galaxies in the Milky Way at My > — 4, 
which indicates that to fully explain the missing satel- 
lites problem, star formation before reionization needs 
to be suppressed, in addition to the — now standard — 
suppression by photo-heating after reionization. 

We have shown by following the linear evolution of 
density fluctuations while taking into account the non- 
linear effect of the supersonic baryon-CDM flow, that 
the stellar-mass function before reionization (z « 11) is 
significantly suppressed at the low-luminosity end, with 
a typical suppression by 50 percent for luminosities today 
of My > —4, going up to 90 percent for regions with a 
high initial supersonic relative velocity. This suppression 
is mainly the result of the reduced accretion of gas onto 
low-mass halos due to the supersonic flow, with a small 
contribution from the overall suppression of the halo 
mass function at 10 5 < M ha i o /M < 10 7 . While these 
calculations use the linear evolution equations for per- 
turbations, cosmological high-resolution hydrodynami- 
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cal simulations of the formation of the first structures 
have shown that the approximations of §$2] and [3] ad- 
equately describe the formation and baryon fraction of 
the fi rst galaxies (jNaoz et al.ll2009l : lO'Learv fc McQuinnl 
l2012f ). potentially even underestimating the suppression 
in the star -forming gas fraction at t he high-mass end 
(« 10 7 M^: iMcQuinn fc Q^Leax^l^oTl ). 

To determine whether the effect of the supersonic 
baryon-CDM flow on the abundance of luminous low- 
mass galaxies persists to the satellite luminosity function 
of a Milky- Way sized halo today, we have run extended 
Press-Schechter simulations of the merger history of a 
1O 12 M halo, taking into account the effect of the super- 
sonic flow on the halo mass function and star formation 
prior to reionization. We found that the effect largely 
remains the same and that the number of satellite galax- 
ies with My > —4 in a Milky- Way sized halo is typically 
suppressed by 50 percent. When then predicting the ob- 
served number of faint satellites of the Milky Way that 
could have been observed by the SDSS — multiplying the 
predicted counts by the fraction of the effective volume 
of the Milky Way observed by the SDSS — we find that 
a typical initial relative velocity of w 30 km s _1 allevi- 
ates most of the discrepancy at the faint end (My > —6) 
between the observed number of satellites and that pre- 
dicted by a model without pre-reionization suppression 
of star formation, without applying any additional mass- 
dependent reduction in star- formation efficiency, e.g., 
due to the radiative feedback from the first stars. There- 
fore, the effect of the supersonic baryon-CDM flow nat- 
urally provides the amount and mass-dependence of the 
pre-reionization suppression of star-formation efficiency 
needed to explain the observed luminosity function of 
Milky- Way satellites. The limited detection efficiency of 
the SDSS for low-luminosity galaxies is such that satel- 
lites at My = — 2 can only be detected out to approx- 
imately 40 kpc, such that a large fraction of the virial 
volume of the Milky Way remains unexplored at these 
low-luminosities. Next-generation surveys such as the 
LSST will be able to observe My = s atellite galaxies 
out t he virial radius of the Milky Way (jTollerud et al.l 
I2QQ8[ ). and the effect of the supersonic baryon-CDM flow 
should lead to a clear suppression in the number of low- 
luminosity satellites. 
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